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Abstract 

In this paper, by making use of the perturbative expansion around 
topological field theory we are trying to understand why the standard 
perturbation theory for General Relativity, which starts with linearized 
gravity does not see gravitational collapse. We start with investigating 
classical equations of motion. For zero Immirzi parameter the ambi- 
guity of the standard perturbative expansion is reproduced. This am- 
biguity is related to the appearance of the linearized diffeomorphism 
symmetry, which becomes unlinked from the original diffeomorphism 
symmetry. Introducing Immirzi parameter makes it possible to restore 
the link between these two symmetries and thus removes the ambi- 
guity, but at the cost of making classical perturbation theory rather 
intractable. Then we argue that the two main sources of complex- 
ity of perturbation theory, infinite number of degrees of freedom and 
non-trivial curvature of the phase space of General Relativity could be 
disentangled when studying quantum amplitudes. As an illustration 
we consider zero order approximation in quantum perturbation the- 
ory. We identify relevant observables, and sketch their quantization. 
We find some indications that this zero order approximation might be 
described by Doubly Special Relativity. 
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1 Introduction 



One of the simplest models describing a collapse of a spacetime is 2 + 1 di- 
mensional gravity coupled to point particles [1] . A particle induces a conical 
singularity with the deficit angle proportional to its mass. If cosmological 
constant is zero spacetime collapses when deficit angle becomes larger than 
2ir. Physically it results in a bound on energy, or in the appearance of an 
invariant energy scale in the relativistic symmetries of particles propaga- 
tion in such spacetime. The latter are known as Doubly Special Relativity 
symmetry [2] (see also [3].) 

The picture of collapsing spacetime becomes more elaborate if we add 
cosmological constant. If A > spacetime has an outer cosmological horizon. 
When we add a massive particle the cosmological horizon shrinks around it, 
and if the mass of the particle exceeds Planckian the horizon shrinks to a 
point. If A < 0, the energy is not bounded from above. However as soon as 
mass of a particle gets larger than mpi, it becomes surrounded by a horizon 
known as the BTZ black hole. 

In four spacetime dimensions the picture of gravitational collapse is even 
more intricate. Classically, any point particles creates a black hole around it 
with its size proportional to its mass. In the presence of positive cosmological 
constant one has also a cosmological horizon which shrinks as the mass 
grows. One can even imagine a situation when the horizon of a black hole 
collides with cosmological horizon. This could be a bound on energy in four 
dimensional GR analogous to the energy bound in 2+1 gravity mentioned 
above. 

It is worth mentioning here that as it is understood at least in 2 + 1 
gravity the phenomena described above are vital for consistent quantization 
of the theory. The energy bound is responsible for the appearance of fixed 
point in renormalization group. 1 

The problem is that the standard perturbative approach which starts 
by linearizing gravity around a fixed background can see none of these phe- 
nomena. If in 2 + 1 dimensions one can hope to solve the theory exactly, 
this seems hardly possible in dimension four. There the complexity of the 
theory is multiplied by the presence of infinite number of degrees of freedom. 
The natural question then arises: is it possible to develop an approximate 
method to solving (quantum) Einstein's equations which would be able to 

1 Sometimes it is said that 2 + 1 gravity owes its renormalizibility to the absence of 
local degrees of freedom. This is not the case however. A simple counterexample is a 
quantization of 2 + 1 gravity by expanding around a fixed background. It knows about 
the absence of propagating modes of gravitational field, but it is still non-renormalizible 
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see gravitational collapse, maybe in some simplified form? The first step 
towards answering this question would be to understand the reason why the 
standard perturbative approach does not manage to do it. 

We try to address this question in the present paper. For this we uti- 
lize the perturbative framework of [4] which uses a certain topological field 
theory (namely BF theory) as a starting point. It requires a non-zero cos- 
mological constant and also admits inclusion of Immirzi parameter. 

We start by studying such perturbative expansion of classical equations 
of motion with zero Immirzi parameter. It appears that it reproduces the 
results of standard perturbation theory which starts with linearization. The 
linearized Einstein's equations emerge in it as certain integrability conditions 
for the first order equations. 

The immediate problem that appears at this stage is the invariance of this 
equations with respect to linearized diffeomorphisms which become unlinked 
from the actual diffeomorphisms and therefore do not correspond to any 
physical symmetry. This leads to ambiguities in the physical results obtained 
from such perturbation theory. 

This may seem strange for an approach starting from BF theory as it 
is invariant with respect to exact diffeomorphisms. However linearized dif- 
feomorphism transformation now appear in a different guise: as a part of 
"translational" symmetry group of BF-theory which is dual to the local 
gauge symmetry. 

A small modification of BF theory, namely introducing Immirzi param- 
eter provides a natural solution to this problem. It connects the "trans- 
lational" symmetry and the local gauge symmetry into a single one, and 
thus, to first approximation, reestablishes the link between linearized and 
actual diffeomorphisms. It has a physical effect that stringy degrees of free- 
dom coming from breaking "translational" symmetry now become coupled 
to particles degrees of freedom coming from breaking local gauge symmetry. 

However the next problem emerges here. The first order integrability 
conditions turn into the form of exact Einstein's equations and as such can- 
not be solved in a generic situations. This makes the classical perturbation 
theory not useful in practice. 

Fortunately the situation is very different in quantum theory where we 
are interested not in calculating the values of the fields everywhere in space, 
but only in scattering amplitudes for finite number of particles. At a fi- 
nite order of perturbation theory only finite number of modes of the field 
contribute to such amplitudes. This allows us to reduce the problem to 
finite-dimensional, while keeping track of some non-linearities of General 
Relativity at the same time. 
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One of the effects of introducing Immirzi parameter in this case is that 
quantum perturbation theory becomes "semi-holographical" . Its Feynman 
diagrams map on a boundary of space and they braid as in 3 dimensions. 
At the same time every particle in a diagram has 3 degrees of freedom as 
in 4 dimensions. In Euclidian theory for example the braiding effect makes 
energy and momenta bounded from above. The bound is controlled by a 
certain combination of Newton constant, cosmological constant, and Immirzi 
parameter. 

In section 2 we review the setup of [4] and [5] for perturbation theory 
for gravity and its coupling to particles. 

In section 3 we solve it to the first order for zero Immirzi parameter and 
identify the problem related to the appearance of an extra gauge symmetry. 

In section 4 we discuss physical implication of this extra gauge symmetry, 
which is decoupling of some of the degrees of freedom linked in the full 
theory, and show how introducing Immirzi parameter restores this coupling. 
We also show that Immirzi parameter mixes mass with Taub-NUT charge, 
which results in the appearance of non-local singularities - gravitational 
Dirac strings. 

In section 5 we give a brief perspective on how to calculate particle scat- 
tering amplitudes in this perturbation theory. We identify the observables 
which contribute to perturbative calculation of the amplitudes and show 
that at any finite order there is only a finite number of them, which makes 
quantum perturbation theory tractable. For non-zero Immirzi parameter 
this observables become non-local due to Dirac strings, and we argue that 
this leads to a certain kind of "dimensional reduction" in quantum ampli- 
tudes. Finally, we take a closer look at the observables which are relevant at 
zero order of perturbation theory. We show that the variables canonically 
conjugate to particle positions form de Sitter space, which is an indication 
of the presence of Doubly Special Relativity effects in the theory. 

In section 6 we discuss some mathematical problems which has to be 
solved to complete the definition of these perturbation theory. 

2 Equations for gravity with particle source 

To set the stage, let us recall, following [4] and [5] the formulation of gravity 
as a constrained topological field theory and it coupling to point particle(s). 
Let us start with the pure gravitational field, which is encoded in a S0(4, 1) 
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(de Sitter) connection 

A M = A IJ T U = Qe/ T a4 + ^ ab T ab ^j (1) 

where index a runs from to 3. In this formula T/j, I, J = 0, . . . 4 are 
generators of Lie algebra so(4, 1). e^ a is the frame field from which the 
metric is constructed and uj fl ab is the spin connection. Connection has a 
natural mass dimension 1 whereas the frame field is dimensionless; this is 
the reason why a length scale t appears in the expression of the components 
of A representing the frame field. This length scale has to be the cosmological 
length I related to the cosmological constant by 

1 _ A 

¥ ~ 3 

To formulate the theory we also need a two form valued in the Lie algebra 
so(4, 1) denoted by B = B^ u ij Tjj dx M A dx u . In terms of this two form and 
connection A the action takes the form 

S = J (Bjj AF"-f e IJKLM Bjj A B KL v M - I B IJ A 5„) (2) 

with vm being a constant vector in fundamental representation of the gauge 
group, forcing the gauge symmetry breaking from SO(4, 1) down to SO(3, 1) 
(see [6].) This action can be shown [4] to be equivalent to the action of Gen- 
eral Relativity with nonzero cosmological constant and a nonzero, dimen- 
sionless Immirzi parameter 7. The initial parameters dimensionless a, (3, £ 
are related to the physical ones as follows 

j_ = A = GA^_ GA 7 

e 2 3' a 3 (i + 7 2 )' 3 (1 + 7 2 ) 1 ' 

Let us now turn to the particle(s) coupling. As it was shown in [5] matter 
can arise in the most natural way by introducing the simplest possible term 
breaking the gauge symmetry of the theory in a localized way. The gauge 
degrees of freedom are then promoted to dynamical degree of freedom, which 
reproduce the dynamics of a relativistic particle coupled to gravity. This idea 
is realized is obtained by choosing a worldline P and a fixed element K of 
the so(4, 1) Lie algebra with generators T IJ , depending on the particle rest 
mass and spin 

K = meT 04 + sT 23 (4) 

so as to have 

S P (A) = - /drTr(KA T (r)) (5) 
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where r parameterizes the world line z M (r) and A r (r) = A jU (z(r)) i M . 

The action for the particle is obtained from (5) by adding gauge degrees 
of freedom that become dynamical at the particle worldline. To this end we 
take A h = h" 1 Ah + h _1 <ih to be the gauge transformation of A. Then the 
particle lagrangian takes the simple form 

L(z,h;A) = -Tr (KA^(r)) S = J dr L(z,h; A) (6) 

which can be rewritten as 

L(z, h; A) = — Tr(JA T ) + L±(z, h) (7) 

where in the first term J is given by 

J^hKh- 1 (8) 

The first term in equation (7) describes the covariant coupling between the 
particle and the A connection of the (constrained) BF theory, while the 
second 

Li(*,h) = -TrOi^hK), (9) 

describes the dynamics of the particle. 

The field equations for gravity coupled to point particle, resulting from 
(2) and (7) in the gauge, in which the particle is at rest at the origin of the 
coordinate system take the form 

F IJ = ae IJKLM B JK v M + pB IJ (10) 

V A B IJ = J IJ dp, 5 P = 5 3 (x)e (11) 

where T>& is covariant derivative of connection A and e is the volume three- 
form on constant time surface. 

In this paper we will be primarily interested in solving these equations in 
the limit a — > 0. It would seem that it suffices to put a = in (10) and solve 
the resulting equations; however there could be contributions from higher 
order equations affecting this limit. Therefore more care is needed and we 
consider expansion of the equations up to the first order. 

Consider the expansion in a 

A = A(°) + «AW + . . . , B = B(°> + aB« + . . . 
In the zero order in a we get 

F /J (AW) = (3B^ IJ (12) 
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V m B^ IJ = J IJ (13) 
In the first order in a we find 

V m AW = e IJKLM B^jK vm + PB^ IJ (14) 

%)5 (i)" + [ A W,B(o)]" = o (15) 

If we now take the covariant differential P A (o) of (14) and compare it with 
(15), we find the first integrability condition 

e IJKLM V MO) (B^ JK v M ) = 0. (16) 

The second integrability condition comes from acting by £> A (o) on eq.(15) 
and reads 

[ B (o),bW*] /j + [A( 1 ),J] /j = 0, (17) 

where 

B(°> IJ = € IJKLM B^jkv m (18) 

Equations (16), (17) affect the zero order contributions to the field B( ) 
and cannot be ignored even in the limit a — > 0. 

Thus the topological limit of gravity, coupled to a point particle is gov- 
erned by four equations (12), (13), (16), (17). In what follows we will drop 
the (°) superscript to simplify the notation. 



3 Zero Immirzi parameter and ambiguities in the 
solution 

In this section we solve the above equations to the first order for (3 = and 
identify the problem arising in such approximation. 

Zero order is described by equations (12), (13). For (3 = those equations 
possess the following gauge symmetry: 

A g- 1 Ag + g^dg, B g- 1 Bg (19) 

and 

B^B + £> A $ (20) 

where g G SO(4, 1), and <E> is a so(4, l)-valued one form. 
The solution in the zeroth order is thus 

A = g-^g B = B p + g-^g, (21) 
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where Bp is a particular solution depending on J. We consider spinless 
particle, so we choose J = g~ 1 K g = m£g _1 T 04 g. 

If we define B p = gB p g _1 the above equation simplifies to 

dB p = K (22) 

and a possible solution is 



Bp = mi T 04 e abc ^ dx b A dx c (23) 



x 



r 



So far g and were arbitrary, but from the first order we have extra 
equations (16), (17). The origin of this equations can also be explained 
the following way. We take the general solution of zeroth order equation 
(21) and plug it into the symmetry breaking term in (2) (the one which is 
proportional to a). Then by varying the resulting expression w.r.t. g we 
obtain equation (17) and by varying it w.r.t. $ we obtain (16). 

However because the whole action is diffeomorphism invariant and be- 
cause the translational part of g is closely related to diffeomorphism trans- 
formations we can start by making an arbitrary choice of g which will simply 
result in a choice of coordinate system. 

For spherical static coordinates in deSitter space we have to choose 



/ cosh | I sirm f yjl-^ sinh \ \ 

J\ — jz sin 9 cos 4> — sin <fi — cos 9 cos <j) —j sin 9 cos 4> 



\l 1 — ji sin 9 sin 4> costp — cos#sin^> — ?sin#sini 



^/1-gcosfl sin# -£cos# 



^ sinh \ j cosh | \j^- jz cosh \ j 

(24) 

It remains now to plug the solution (23) to first integrability condition 
(16) and solve it for For [3 = V~ k in (16) becomes ordinary differential d, 
and the equation reminds the linearized Einstein equation, but with unusual 
type of source. In this case the solution for $ reads 



r 2 \ x b 



* = -e\i- 7i U ijbV3 dtlV 



1 - ^ e iab ^ dx a (sinh j T Ai + cosh t T ^ (25) 
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Notice that to this point the presence of non-zero cosmological constant was 
required. For A — > one has I — > oo and the expression (25) for $ diverges. 

Having solved the condition (16) one can turn to finding from (14). 
The later equation has a regular limit as I — > oo and we can take it at this 
point. We will write down explicitly only the metrical part of it, i.e. the 
part proportional to T 4a 

= ~T m dt + ^T ii dx i + {...)T ab + V m V, (26) 

where ^ is an arbitrary 0-form field. This expression contains the standard 
Newtonian potential term, but there is also some ambiguity in it. Different 
choices of VP lead to different metric that cannot be related by coordinate 
transformation. This is a source of ambiguity in physical results. 

Similar problem can be observed in linearized gravity as well. Linearized 
Einstein's equations 

V m V m h ab - 2V m V ( A )m + V a V b ( 5 mn /i mn ) = (27) 

are invariant with respect to exact diffeomorphisms when one simultane- 
ously transforms the background field <?q 6 —> ff™ n ff— §f- an d the fluctu- 
ation h a b — > ^mn^-^r- At the same time they are invariant with re- 
spect to linearized diffeomorphisms when one transforms fluctuations only 
h a b — ► h a b + V( £j,), where £ a is an increment of a coordinate x a , keeping 
the background field fixed. Linearized gravity doesn't "remember" that the 
two above symmetries have common origin. Therefore the gauge-fixing for 
both symmetries can be done independently. 

To illustrate one of the consequences of the above ambiguity let us recall 
one curious fact that the Schwarzschild solution of the Einstein's equations 
in Kerr-Schild parameterization is also a solution of linearized Einstein's 
equations (27) [7]. It is related to the usual Schwarzschild coordinates by 
transformation T = t — 2mln(r/2m — 1), and the metric in this parameter- 
ization looks like 

ds 2 = -(l - 2 -^)dT 2 - ^dTdr + (l + ^)dr 2 + r 2 dQ. (28) 

It is a solution of equations (27), and it looks quite similar to the standard 
Newtonian metric 

ds 2 = -(l - 2 ^)dT 2 + (l + ^)dr 2 + r 2 dn, (29) 
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which is also a solution of (27). The two however cannot be related by 
a coordinate transformation, and they are physically different with differ- 
ent Penrose diagrams. They can be related by linearized diffeomorphism 
transformation with £o = — 4mln(r/4m) and £j = 0. 

Clearly the arbitrariness in the choice of * in (26) leads to exactly the 
same ambiguity as above. By adjusting we can have either standard 
newtonian potential or a black hole. This is a problem because perturbation 
theory cannot distinguish between the two. 

However we also have to remember about the second integrability con- 
dition (17). By direct calculation one can check that away from particle 
sources it is automatically satisfied. This was to be expected as a result of 
diffeomorphism invariance. At the points of the particle sources this condi- 
tion reduces to 

[A« + $*,J] /J = 0, (30) 

where 

$* IJ = e IJKLM <S>j K v M (31) 

As it was shown in [5] equation (30) is directly related to geodesic equa- 
tion for the test particles. 

One can see that in equation (30) the problem of invariance of equation 
(14) with respect to linearized gauge transformation A* 1 ) — > A* 1 ) + T>a^ is 
combined with another ambiguity which is the invariance of equation (16) 
for ^ with respect to similar transformation 

<D - $ + V AX , (32) 

where x is an arbitrary scalar field. But the later symmetry is a subgroup 
of the symmetry of .BF-theory, namely (20). Thus, the symmetry that 
causes the problem in the linearized gravity has its counterpart in BF theory. 
Therefore, the perturbative approach starting with BF theory considered in 
this section will suffer from the same problem as the standard perturba- 
tion theory that starts from linearized gravity. Similar conclusion has been 
reached also in [8] in the case of an analogous treatment of Yang- Mills theory. 

However the considerations of this section provide a clear hint on how 
this problem is to be cured within perturbative framework. This is the 
subject of the next section. 
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4 Non-zero Immirzi parameter and restoration of 
unique gauge symmetry 

Let us first discuss the physical meaning of the gauge symmetries that ap- 
pear in BF theory. Initially the gauge parameters physically irrelevant as 
such, of course. However, after introducing matter by applying the symme- 
try breaking perturbation some of the gauge parameters turn into physical 
degrees of freedom. The easiest way to look at them is to introduce localized 
symmetry breaking terms, which results in the appearance of finite number 
of degrees of freedom, similar to those we used introducing particles in sec- 
tion 1. Then one can ask: what is the physical meaning of the degrees of 
freedom resulting from such symmetry breaking. 

Consider BF theory coupled to a most general type of localized sources, 
which includes particles considered above as well as strings considered in [9] . 

S = J Bij A F IJ + S P {k) + S S (B), (33) 
where Sp(A) is the same as in (5), and the string action reads 

S S (B) = J drdaTi(lB Tf7 ), (34) 

where r and a are variables parameterizing the string worldsheet, B ra = 
B[mv tIt" i and L is a charge analogous to (4) but so far with different mass 
and spin. 

Equations of motion for particles were described in section 2. To obtain 
equations of motion for stringy source one has to perform a gauge transfor- 
mation B — > h~ 1 (B + d^)h and substitute it back to the string action (34). 
After integration by part the string action will look like 

S s = f dTdaTi(PB ra ) + / drTr($P) - / dTdaTr(<Z>d A P), (35) 
JW JdW JW 

where P = hl_h -1 is the string momentum density. The equations following 
from action (35) have to be supplied with boundary conditions at the string 
endpoints. To allow for non-vanishing string momentum density the natural 
condition is 

<*> mw = V AX . (36) 
aW 

Thus, from the third term in (35) we get conservation of string momentum 
density along the whole worldsheet, 

d A P = 0, (37) 
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and from the second term in (35) we get the same condition, but restricted 
to string endpoints. The variational principle thus become consistent for 
arbitrary string momentum P. From what the second term in (35) with 
boundary condition (36) turns into one can deduce that the gauge parameter 
X is canonically conjugate to the string momentum P at its endpoints. Notice 
that the gauge parameter \ can be turned into a physical degree of freedom 
only at the endpoint of the string. Thus, it becomes somewhat analogous 
to Chan-Patton degrees of freedom in string theory. 

So far the particle positions h and momenta J are completely indepen- 
dent from string momenta P and their canonical conjugates \. This is a 
consequence of a split between exact gauge transformations encoded in h 
and linearized gauge transformations encoded in \- Below we consider how 
the link between the two can be restored already at zero order of perturba- 
tion theory. 

Consider the constraints of the theory described by action (33). These 
are equations (12), (13) reduced on the spacial slice £ with sources on the 
r.h.s. (so far we are considering (3 = case). 

F IJ (A) = P IJ 5 S (38) 

where 5s = Jwnz, da5 3 (x — x(a)), and 

V A B IJ = J IJ S P , (39) 

where J IJ is a particle source described in section 2. Constraint (39) gener- 
ates gauge transformations (19), and constraint (38) - transformations (20). 
The later contains a subgroup (32), which is directly linearized diffeomor- 
phism transformations that were detected in the previous section as a source 
of the problem. 

The obvious approach to solving this problem is to make the theory 
remember that the two symmetries have common origin, i.e. that the con- 
straints (38) and (39) generating them are not actually independent. This 
is what naturally happens when we introduce a (3 term in (38) 

F IJ (A) + @B IJ = P IJ 5 S (40) 

(39) becomes a consequence of (40) due to Bianchi identity, provided that 
particle sources are located at the endpoints of the strings. The string be- 
comes something like the Dirac string for "magnetic" charge J. The role of 
the Dirac strings in this context we will discuss below in detail. Positions of 
the particles and positions of the strings are thus no longer independent. At 
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the same time the gauge transformations (20) get deformed and also start 
affecting connection. The constraint (40) generates the following transfor- 
mations 

B-> B + 2? A $ + A^A + /3$, (41) 

and the analog of the subgroup (32) is now 

A + / 9$^g- 1 dg + g- 1 (A + /3$)g, (42) 

which is the same as the gauge transformation (19). The link between 
the original gauge transformation and the linearized gauge transformation 
is thus restored. One of the resulting effects is that particles degrees of 
freedom and Chan-Patton degrees of freedom at the strings endpoints have 
merged into one. 

Let us now take a closer look at equations (12), (13), and (14) for (3 / 
0. The first, rather obvious thing to notice is that, as a result of Bianchi 
identity, eqs. (12) and (13) cannot be solved for non-trivial source if the 
connection is nonsingular. This problem can be solved by recalling the 
Dirac magnetic monopole solution (see, for example, [10]): we assume that 
our connection contains the part proportional to Dirac monopole connection 

A = g- l dg + (3g~ l {A D + <5>)g, A D ~a D = (l-cos0)# (43) 

where g 6 SO(4, 1), and $ is a so(4, l)-valued one form. We will use the 
ansatz (43) below. In fact the Dirac string and the string described in the 
beginning of this section are the same. 

Consider now eq. (16), which is an integrability condition for the first 
order equations. Computing covariant differential, and noticing that since 
vm has only one non-zero component i>4 = 1 we find 

Pa"m = A M 4 = ej, i = 0, ...4 

where is the tetrad one form. Now it follows from eqs. (16) and (12), (13) 
that 

^e IJKL *F JK Ae l + e IJKLM v M J KL = 

It can be shown that for one particle the second term in the expression above 
vanishes, and thus we are left with 

e IJKLi FjKAe . = Q (44) 
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Decomposing F into torsion = Tj and curvature Fij = Rij two-forms one 
finds that any connection being a solution of our problem must be torsion- 
free Tj = and satisfy the vacuum Einstein-de Sitter equations 2 

e ijkl R jk Ae t = (45) 

Having obtained the exact Einstein's equations already at the first or- 
der of perturbation theory, one can be sure that non-perturbative effects, 
including gravitational collapse cannot be missed. But this is not a solu- 
tion of any problem, because doing such perturbation theory is the same 
as solving General Relativity non-perturbatively. Thus, in classical theory, 
the present perturbative approach does not lead to any simplification. How- 
ever it can bring about some new intuition in quantum theory which is the 
subject of the next section. 

Before finishing this section let us point out some another consequence 
of introducing Immirzi parameter. Let us solve eq.(45) for the ansatz (43) 
to the leading order in f3, which is the same as linear approximation. In this 
approximation the solution for $ will be (25). One can check that g _1 <I>g 
does not have terms proportional to T 4c so that the contribution of it to the 
tetrad one-form is zero. The only contribution to the tetrad (in the leading 
order in (3) is therefore (cf. (43) 

A i4 = e i= (g- 1 dg + Ja D ) (46) 
One can then easily find what the metric g - 




9 = -{ !-^J {dt + N{l-cos6)d(t>) 2 +[l-^J dr 2 +r 2 (d9 2 +sm 2 9 d(p 2 ). 

(47) 

This is de Sitter-Taub-NUT metric linearized with the Taub-NUT charge N 
and with no mass. If one goes beyond linear approximation one obtains the 
exact Taub-NUT solution. This means that the Immirzi parameter mixes 
mass with Taub-NUT charge. The presence of a non-local (Dirac string) 
singularity has a profound consequences in quantum theory. This will be 
discussed in the next section. 



5 Quantum amplitudes versus classical solutions 

The two main sources of difficulties with four dimensional General Relativity 
are infinite number of degrees of freedom and non-linearity of its equations. 

2 That is, Einstein equations with a positive cosmological constant. 
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One of the consequences of the latter is a non-trivial curvature of the phase 
space of the theory. One of the lessons from 2+1 gravity coupled to matter is 
that momentum part of its phase space is a group manifold rather than linear 
space and the information about its curvature is essential for consistent 
quantization of the theory. 

In standard perturbation theory one starts by choosing a background 
- a certain point in the phase space of the theory. One then proceeds by 
linearizing the theory around this background, which in particular implies 
that one works in a tangent space of the phase space of the theory instead 
of the actual phase space. This space remains linear even after removing 
gauge degrees of freedom by taking its quotient by linearized diffeomorphism 
transformations. 

Another property of quantum perturbation theory is that when we are 
interested in scattering of a finite number of particles in a finite order of per- 
turbation theory we have to calculate Feynman diagrams. The expressions 
for Feynman diagrams are finite dimensional integrals instead of infinite di- 
mensional path integrals. Thus quantum theory allows us to select a set of 
questions which could be answered by taking into account only finite number 
of modes of quantum field. This option is not available in classical pertur- 
bation theory where we have to solve infinite-dimensional field equations in 
every order. This is a big advantage of quantum perturbation theory over 
the classical one. 

Now one can notice that the ability of quantum perturbation theory to 
reduce the problem of calculating scattering amplitudes to finite dimensional 
one, and its inability to see the curvature of the phase space of the theory 
are not logically linked which each other. Therefore one can ask if there 
is a quantum perturbation theory which allows to reduce the problem to a 
certain finite-dimensional subspace of the phase space and if that subspace 
is curved the perturbation theory still can see it. 

Below we show that the quantum analog of the perturbation theory 
studied here does precisely this job, and it allows us to see some non-trivial 
effect already in zeroth order. 

We start by explaining how the perturbative expressions for quantum 
amplitudes get reduced to finite-dimensional integrals. 

The general perturbative expression for the partition function coupled 
to arbitrary finite number of particles (where we neglect all the interactions 
except gravitational) after integrating out S-field looks like 

Z({x Pi },{x Pf }) = [ VAJ2 ( -^j-( I v A e ABCDE F BC (x)AF DE (x)) n 
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xexp[i f F IJ AFjj + Y^Spix^Xp )], (48) 

JM 

Where S p are actions for particles of the type (5) , and x Pi and Xpj. 9)1*6 initial 
and final points of p-th. particle. They are presumably located on initial and 
final slice of M. The expression (48) can be rewritten in terms of n-point 
functions as 

Z({x Pi },{x Pf }) = ^—y- ( dx 1 ...dx n W(x 1 ,...,x n ,{x Pi },{x Pf }), (49) 
where 

W{ Xl ,...,x n ,{x Pi },{x Pf })= f VAH(v A e ABCDE F BC (x n ) A F DE (x n )) 

exp[4 / F IJ AFu + Y, S P (x Pi ,x Pf )} (50) 

P J M 

Now one has to recall one of the key properties of n-point functions of 
diffeomorphism invariant theories, i.e. when the action and the measure are 
diffeomorphism invariant, (see e.g. [11] and references therein) that they 
do not actually depend on x n , as coordinates do not carry any physical 
information. The only dependence may come from possible coincidences of 
x n , but it yields a measure zero contribution. Of course n-point functions 
depend on relative positions of the points x n , but the latter are encoded not 
in the coordinates but in the fields. 

This means that the integrands in (49) are constants and therefore the 
integrals can simply be dropped, i.e. absorbed into renormalization of a. 
As a consequence, the n-order contribution to the partition function can 
be evaluated as a path integral of topological field theory with n localized 
symmetry breaking insertions. As it was discussed e.g. in [5]: breaking a 
symmetry of topological field theory in a localized way releases only a finite 
number of degrees of freedom. This is the central observation which allows 
us to reduce the procedure of perturbative calculation of quantum scattering 
amplitudes to a finite dimensional problem. Notice that this would not be 
possible in non-perturbative calculation as in such case the integral with a 
symmetry breaking term would stay in the exponent and as a result could 
not be replaced with a local symmetry breaking insertion. 

Now let us try to compose a list of variables which become relevant 
for the calculation of the n-point function (50). To simplify the job let us 
first integrate out the particles degrees of freedom, g p , which is the gauge 
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parameter evaluated at the location of the particle. It can be done by 
replacing particle terms in the action with particle propagator. 

G p {g(x Pi ),g(x Pf ), A p ) = J Vgexp ( J <ftTr(A t J) + S p {g)) (51) 

We will not dwell on calculating this expression here. However, as it was 
stressed in [5], the propagator (51) is a matrix element of the open Wilson 
line operator whose representation is defined by the mass and the spin of the 
particle. The states, between which the matrix element has to be taken are 
SO(4, 1) analogs of Wigner functions, with corresponding representation 
indices. Alternatively, one can insert Wilson lines directly. This would 
correspond to working in momentum representation, as the matrix indices 
of the Wilson line are momenta of the particles. Eq. (50) now reads: 

W(xi, ...,x n ,{x Pi },{x Pf }) = / VA Y[ (G p (g(x Pi ),g(x Pf ))) 

J p 

H(v A e ABCDE F BC (x n )F DE (x n ))eMi [ F IJ A F u ], (52) 

n Jm 

The dependence of this expression on {x n }, {x Pi }, and {x Pf } should be 
understood as the dependence on the connection field evaluated at these 
points. 

Now if instead of connection field A we use its holonomies connecting dif- 
ferent point in space, it is clear that due to 50(4, 1) gauge symmetry of the 
action in the exponent the integrand in (52) cannot depend on holonomies 
ending anywhere else but at the points x Pi , x Pf , and x n . For convenience 
however let us choose one common reference point x$ putting it for exam- 
ple somewhere on the boundary of initial slice, xo G dM n Sj. For basis 
holonomies we choose g n connecting x n to xq, g Pi connecting x Pi to xq, and 
g Pf connecting x Pf to Xq. Any holonomy connecting any pair of points can 
be expressed as a combination of those from the above list. If we took 
f3 = limit, the free action would force the connection to be flat, and the 
holonomies would become independent of pathes along which they are taken. 
That means that these holonomies would comprise the full list of variables 
relevant for calculating the amplitude (52). 

The situation is very different when we consider (3 ^ 0. First of all as it 
was discussed in [5] the BF theory with nonzero (3 cannot be consistently 
coupled to a local symmetry breaking source without introducing introduc- 
ing an extended connection singularity - the Dirac string. 
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Generally, if in classical theory the equations of motion are inconsistent 
because of a charge non-conservation, in quantum theory in analogous sit- 
uation one obtains a vanishing partition function. The situation is similar 
here. For example, if a particle located away from the boundary, and the 
connection is non-singular, in the partition function 

Z = J VAexp (i J F IJ A F u + iTr(AK)5 3 (x - z(t))) (53) 

The first term in the exponent will be independent of A(z(t)), and the con- 
tribution from integrating over it will be of the form S(K). This means that 
if K 7^ the partition function is zero, and the only way to get non-zero 
contribution to the partition function is to either introduce the Dirac string, 
or put the particle on a spatial boundary. Below we will argue that the two 
are equivalent. 

For similar reason the local symmetry breaking insertion, say at the point 
x, also has to be connected by Dirac string to the boundary. The easiest 
way to see this is to extract an integral over dg(x) from the measure VA, 
where g{x) is some holonomy connecting x to the reference point along some 
path. If DA is the Ashtekar-Lewandowski measure [12], it has to contain 
dg{x). For simplicity consider 1-point function. The argument however will 
hold for n-point functions for non-coincident x n as well. We will extract 
explicitly g{x) variable, after which the 1-point function can be written as 

W(x)= J T>Adg(x)Tr(g(x)~ 1 vg(x)F(x) A F(x)) exp (i J F IJ A Fijj , (54) 

where v is the matrix representation of the vector v A , in an appropriate 
representation. Now one can notice that g(x) is a functional of connection 
in the bulk, while the exponent term in the absence of Dirac string does 
not depend on the bulk connection. This means that the only dependence 
of the integrand in (54) on g is in pre-exponential term. The integral over 
dg{x) can be therefore explicitly calculated. It equals zero because it is 
an averaging of a vector over all possible directions. The only way to get 
a nonzero contribution is to introduce the Dirac string as in the case of 
particle. 

The presence of the Dirac strings makes the holonomies path-dependent: 
the holonomy passing to the left from the string will be different from the 
one passing to the right. To keep track of this information we should to the 
list of all holonomies relevant for the (3 = case described above add those 
around each string. A holonomy around the string ending at x n , x Pi , or x Pf 
with the initial point connected to the tip of the string will be called g*., 
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and g* respectively. The natural choice of pathes for holonomes g n , g Pi and 
g Pf will be then the following. g n for example will start at x n which is the 
tip of one Dirac string, then follow this string towards the boundary, and 
then go along the boundary to the reference point xq. 

One can easily see that the above picture is in fact three dimensional 
instead of four. A Dirac string (more precisely a thin tube cut around it) can 
be considered as an extension of the boundary of space. The tip of each Dirac 
string can be considered as a puncture on the boundary. The holonomies 
along the strings can be then considered as holonomies taken along the 
boundary and connecting punctures. And the holonomies around strings 
just become holonomies around punctures. This dimensional reduction is, 
of course, just a reflection of the fact that the action in the exponent in (52) 
is the action of three dimensional theory, namely the Chern-Simons theory. 

In the present paper we will be interested only in the limit a — > 0. 
Quantum theory allows us to consistently take that limit just by dropping all 
the terms proportional to a in the expansion (49), no care about integrability 
conditions as those encountered in classical theory is needed. 

The relevant n-point (actually 0-point) function will simply read: 

W({x Pi },{x Pf }) = f VA]J(W p (x Pi ,x Pf ))e W [i I Y CS {A)\, (55) 

J p J 8M 

Following the discussion above in this expression we have taken into account 
that all the points x Pi and x Pf and all the Wilson lines connecting them can 
be mapped on the three dimensional boundary of space. 

Instead of doing path integral we will try to deduce physics by making 
canonical analysis of the theory encoded in (55). For this it is enough to 
consider the variable residing on the initial slice of dM, ie those connecting 
x Pi . Let x Po be a reference point, and x Pl is the particle of interest. The 
position of the particle with respect to the reference point will be described 
by an SO (A, 1) group element - the holonomy goi connecting x Po and x Pl . 
Notice that although the expression for the amplitudes is described in terms 
of three dimensional theory the particle still has 3 translational degrees of 
freedom as does four dimensional particle. 

Now one can ask what is the momentum of the particle, i.e. what is 
the variable canonically conjugate to translational part of go\. The Poisson 
structure of Chern-Simons theory 

{Af B ,Af D } = (3e ij 6t AC 6 B ^, (56) 
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tells us that it has to be the translational part of the holonomy g\ which 
is taken around the particle x Pl which starts and ends at x Po . g\ is also 
an element of SO{A, 1) group, so the momentum space of the particle is de 
Sitter space. This strongly indicates that a theory describing relativistic 
symmetries of such a system of particles, in the flat space, A — > limit will 
be of DSR type, like in the 2+1 dimensional case. 

This expectation is also motivated by strong results on Hamiltonian 
quantization of Chern-Simon theories with arbitrary gauge group on Rie- 
mann surfaces [13]. On the basis of these results one expects that after 
quantization the symmetry group SO(4, 1) will be replaced by its quantum 
counter part U q {SO{A, 1)) (Note however that Fock-Rosly theorem applies, 
strictly speaking, to compact gauge group and its extension to the non- 
compact case is highly nontrivial, see [14] for the analysis of SL(2, C) case.) 
Then the limit of cosmological constant going to zero will presumably lead 
to emergence of «;-Poincare algebra [15], which is central for DSR program. 

To finish this section let us say a few words about euclidian version of 
this model. For SO (5) gauge group the momentum space is a sphere, both 
energy and momenta are bounded, and the theory has a natural UV cutoff. 
Spacetime collapses at the energy higher than y/A//3. It is not yet clear 
what is the precise mechanism of this collapse. There are two possibilities 
in principle. Either all spacetime disappears at one instant, as it is the case 
for 2 + 1 gravity with zero cosmological constant. Or it starts collapsing 
in some region around a particle, and the region will be growing as the 
mass of the particle is growing. In the later case there may be some effects 
not suppressed by the smallness of the cosmological constant. This will be 
studied in the subsequent papers. 

6 Discussion 

To this point we do not have explicit expressions for higher order corrections 
to scattering amplitudes. In this section we will briefly mention the problems 
one has to face in calculating them. 

Let us try to reexpress (52) in terms of finite dimensional integrals over 
holonomies instead of path integral over connections. First one can replace 
the action in the exponent in (52) by Chern-Simons action. Then we can take 
a locally flat connection connection A ({g}) whose moduli space is defined 
by holonomies g n , g Pi , g Pf and <?*, g*., g* f (we will denote the whole set by 
{g} ), and plug it into the Chern -Simons action. We will have a functional 
depending on only finite number of parameters. 
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I Y cs (A({g})) = S CS ({ 9 }) (57) 

A more complicated treatment is needed for preexponent terms. The lo- 
cal curvature has somehow to be read of from the set of holonomies. Clearly 
the two relevant holonomies at a point x n are g n and g*. They are dual 
to each other and therefore the should give rise to two dual elements of 
curvature entering a local Fbc(x u ) A Ff)E(x n ) insertion in the action. The 
difficulty here is that holonomies are group elements and curvature is an 
algebra element. Therefore we need a certain map J- IJ (g) from from Lie 
group to Lie algebra. Having such a map we can substitute in our path 
integral 

FbcM A F DE (x n ) - T BC (g n )T DE (g* n ) (58) 

to obtain 

W(x 1 ,...,x n ,{x Pi },{x Pf }) = / Y[{g}dgY[{G p (g~^g Pf )) 

J p 

II (v A e ABCDE T BC {g n )T DE (g* n )) expfe(M)], (59) 

n 

This is a "Feynman diagram" expression for an n-point function in this 
perturbation theory. The only missing ingredient is the exact form of the 
map T . The only natural requirement on this map is gauge invariance 

T 1J (h- x gh) = h^^ht (60) 

This however does not specify the map uniquely and it has to be derived 
from other requirements. Our proposal would be to derive the exact form 
of T by comparing an analog of the expression (59) but without symmetry 
breaking with derivatives of Crane -Yetter invariants with respect to its 
coupling constant. 
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